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. We propose an operational definition of the entropy of cosmological perturbations based 

on a truncation of the hierarchy of Green functions. The value of the entropy is unambiguous 
despite gauge invariance and the renormalization procedure. At the first level of truncation, 
' the reduced density matrices are Gaussian and the entropy is the only intrinsic quantity. In 

this case, the quantum-to-classical transition concerns the entanglement of modes of opposite 
wave- vectors, and the threshold of classicality is that of separability. The relations to other 
criteria of classicality are established. We explain why, during inflation, most of these criteria 
are not intrinsic. We complete our analysis by showing that all reduced density matrices can 
be written as statistical mixtures of minimal states, the squeezed properties of which are less 
constrained as the entropy increases. Pointer states therefore appear not to be relevant to 
the discussion. The entropy is calculated for various models in paper II. 



> ■ I. INTRODUCTION 

00 

l/-^ . The standard predictions of inflation are in noteworthy agreement with the recent 

CD \ observations [H . Yet several aspects of inflation remain poorly understood. Among these 

lO ' is the so-called quantum-to-classical transition of cosmological perturbations. 

. In the inflationary scenario, primordial power spectra of gravitational waves and scalar 

\ perturbations result from the parametric amplification of vacuum fluctuations which be- 

gins once the modes exit the horizon. In the course of this unitary evolution, modes of 
. opposite wave- vector become more and more entangled. However, the primordial power 

^ ! spectra are impervious to this entanglement [1,0] because, for this expectation value, the 

\ relative contribution of the quantum correlations is inversely proportional to the amplifi- 

cation factor. Therefore, as far as (today) observational data are concerned, one can safely 
use a stochastic ensemble of growing modes in place of the pure entangled state predicted 
by the quantum treatment. 

Yet, this distinction matters for other observables. In particular, the calculation of the 
entropy and that of backreaction effects (like any kind of radiative corrections) must be 
addressed in the quantum settings. Although these aspects are related, in these notes we 
specifically consider the question of decoherence. We split the presentation in two papers, 
called I and II. In Paper II Q], we calculate of the entropy in various models. In the 
present paper we focus on the formulation of the question. There is no way around this 
first step. Even though we have some notions of general features of decoherence, there 
is no universal description of it. It occurs in a variety of ways, depending on whether 
the system is chaotic [H], fermionic 0], a two-level system 0], or dominated by tunneling 
events [1]. Decoherence must be analyzed case by case. 
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Let us therefore present the specific difficulties regarding cosmological perturbations. 
The first is of a technical nature, the second concerns the consistency of the hypothesis. 
Cosmological perturbations are described by weakly interacting quantum fields propagat- 
ing on a fixed background geometry. One therefore stumbles upon the infinities inherent 
to any interacting QFT. We argue that decoherence can only be properly formulated in 
terms of expectation values of renormalized quantities. That is, reduced density matrices 
are properly defined by a self-consistent truncation of the hierarchy of Green functions, 
rather than by solutions of a master equation. The second difficulty concerns the nature 
of the environment. It will be investigated in more details in Indeed, since one is 
describing adiabatic perturbations, one cannot introduce "external" dynamical degrees of 
freedom that will act as an environment. (If one does so, as in @, [l^l) decoherence obtains 
by construction. But it is based on a mechanism which might not be relevant in inflation, 
and it occurs at a rate which is unknown, both because the coupling strength and the 
statistical properties of the extra degrees of freedom are ad hoc.) The coarse graining 
should be phrased in terms of the properties of the system itself. The Green functions 
offer this possibility. 

In these notes, we truncate the hierarchy of Green functions at the first non trivial level: 
we retain the two-point correlation function, and set to zero all higher order connected 
correlation functions. The reduced density matrices so defined are Gaussian and Homo- 
geneous Density Matrices (GHDM). These factorize in two- mode sectors characterized by 
opposite wave vectors (k, — k). The class of reduced density matrices being defined, we 
can turn to the calculation of the entropy and to the analysis of the quantum-to-classical 
transition. For the entropy, we argue that its value is unambiguous, despite infinities of 
QFT and gauge invariance. We also prove that, during inflation, the entropy is the only 
intrinsic property of GHDM. Indeed, the values of the other quantities (e.g. the number of 
quanta) depend on the choice of the canonical variables which are not univocally defined 
because the frequency is not constant. We calculate the entropy for two classes of models 
in II. 

To address the question of the quantum-to-classical transition, we need to classify the 
reduced density matrices into 'classical' and 'quantum' states. The quantum properties of 
GHDM are linked to the entanglement of modes of opposite wave vectors. The quantum- 
to-classical transition thus occurs when this entanglement is lost (which happens when 
the state is neither pure nor thermalized, but at a sharp frontier in between). This gives 
an operational definition of classical states and of the time of decoherence for GHDM. 

To apply this definition, one must first select creation and destruction operators. The 
latter are well defined only if the Hamiltonian is time independent. This criterion, as well 
as the other criteria discussed in Sec. IVIl are therefore unambiguous during the radiation 
dominated era [3], but not during inflation. Yet, within a given representation, separability 
is a meaningful concept which yields to a physical picture of entanglement. 

This physical interpretation is developed in Sec. IVIl and IVII Al In particular, in Sec. 
IVII A| we link decoherence to the possibility of writting GHDM as statistical mixtures 
of minimal states (displaced squeezed states), the squeezing of which is less and less 
constrained as the entropy increases. In Sec. IVIIBl we argue that pointer states are 
not relevant to describe the decoherence of primordial fluctuations. 

We have organized the paper as follows. Sec. [H] recapitulates the properties of the 
state of the linear perturbations. The coarse graining is defined in Sec. IIIH where the 
properties of the GHDM are also summarized. Classicality is then defined as separability, 
as explained in Sec. |Vl Sections IIIII and IVj aim also at clarifying our previous work 3. In 
Sec. IIV| we identify the von Neumann entropy as the only intrinsic statistical property of 
these states. This establishes that it is unambiguous despite the redundancy of Einstein's 
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equations and the ambiguities from the perturbative renormahzation of Green functions. 
In Sec. IVIl the definition of separable states is compared to three ahernative equivalent 
criteria of classicality, and to one non-equivalent criterion. Our concluding remarks in Sec. 
IVIIBI concern the irrelevance of pointer states to the question of the quantum-to-classical 
transition. 



II. THE STATE OF LINEAR PERTURBATIONS 
A. Settings 

In models with one inflaton field, the dynamics of linearized (scalar and tensor) pertur- 
bations is similar to the evolution of a massless scalar field (p in the background space-time 



111 ]. The latter is a Friedman- Robertson- Walker space-time with fiat spatial section. The 



line element is 

ds"^ = a'^irj) [-drf + Sijdx'dx^] , (1) 

The Hubble parameter H = d\na/dt is slowly evolving. This variation is governed by e, 
the first slow-roll parameter, e = — din H/ din a. We consider scalar perturbations, the 
treatment of tensor perturbation proceeds along similar lines. In this case, ip designates 
the Mukhanov-Sasaki variable which is a linear combination of the infiaton perturbations 
of the gravitational potential. Its Hamiltonian is unique up to a boundary term which 
corresponds to a particular choice of canonical variables. We will come back to this 
important point in Sec. III.C. If we choose the conjugate momentum of 99 to be vr = drj(p, 
the Hamiltonian describing the evolution of linear perturbations is 



'd-'kHk^-i, 




(2) 

where (/Jk (= V'lk) Fourier transform of the field amplitude, and vTk its conjugate 

momentum. The time-dependent function z relates (p to the scalar primordial curvature 



perturbation (defined on hypersurfaces orthogonal to the comoving worldlines jlll ]) 



(^(t,x)=z(t)C(t,x), z{t) = ^^. (3) 

In single field infiation, the sound velocity is = 1. 

We quantize each two-mode system (yjk^'^-k) in the Schrodinger Picture (SP) which 
is best adapted to describe the Gaussian states considered in the following Sections. The 
field amplitude is decomposed at a given time r/in in terms of time-independent creation 
and annihilation operators 

pM = (air + alt) . Mmn) = -^^[\ (air - a^i) , (4a) 



where a™ and a*"'^ verify the commutation relations 

5(3) (k - k') . (5) 



in in\ 
"k ' "k' 
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The free vacuum will be taken to be the Bunch-Davis (BD) vacuum, defined as the state 
which minimizes the Hamiltonian ([2|) in the asymptotic past 

ajf|Oin)=0, for r/in^-oo. (6) 

Alternately, in the Heisenberg Picture (HP), it corresponds to the state with only positive 
frequencies in the asymptotic past. In the decomposition of the field amplitude ^Pkii]) = 
aj^(/9^"(ry) + a^^ip^^*{r]), the mode functions ipJP' verify 

{id^ - k) ^flkr,^.^ = . (7) 

In a de Sitter background, a = —1/Hr], they are given by 



As is well known, there are no linear scalar perturbations in a purely de Sitter background. 
The solution ^ merely serves the purpose to write the modes in a closed form. In the slow- 
roll approximation, the scalar perturbation spectra in the long wavelength limit (krj <^ 1) 
can be inferred from the above solutions by the substitution H Hf^/y/eil where the 
background quantities are evaluated at horizon crossing, i.e. k = a^Hk. 

The free evolution corresponding to the Hamiltonian ([2|) preserves the Gaussianity and 
the purity of the initial state ([6]). In addition, each two mode sector can be analyzed 
independently since the state and the Hamiltonian split into a tensor product and a sum 
respectively, 



l^in(r?)) = 09 l^k,-k(??)) , 
Ik's 

l^k,-k(»?)) = Texp {- y d7?'i7k,-k(?/)) |Oin,(k,-k)) , (9) 
where the tensor product is over half of the wave vectors since |^'k,-k) is a two-mode state. 



B. The covariance matrix 

Because the pairs of modes are statistically independent in the state Q, we consider 
only one such pair and drop the index (k, — k). To use only one formalism throughout the 
paper, from now on we adopt the density matrix notation 

p{^) = \^>){^. (10) 

Since pij]) is Gaussian, its statistical properties are summarized in the one and two- 
times correlation functions. The former vanish identically by statistical homogeneity. 
Many of the second moments vanish as well, also because of statistical homogeneity, namely 
Tr (p<^^j^) = Tr (/Ovr^.j^) = Tr (p^^k'^k) = Tr (p93~k^-k) = 0. Finally, due to the relation 

= (/7_k and similarly for TTk, the remaning covariances can be conveniently condensed 
into a 2 X 2 matrix (instead of 4 x 4) that we shall also call the covariance matrix 
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where the P's are functions of rj and k. In the inflationary phase (approximated by a de 
Sitter evolution), they are given by 

= liW^, ^-k}) = 2Re [y^rdr,y^f] = , (12c) 

where { , } is the anticommutator and we used the notation x = krj. 

We introduce an additional representation of the state ([9|) which clearly displays the 
entanglement between modes of opposite wave vectors: 

n^Tr(^p{r^)a'^^a^) , (13a) 
c^Tt{pi^)a^a'\) , (13b) 

n is real while c is complex. The utility of this represention stems from the fact that n is 
simply related to the power spectrum while |c| measures the strength of the correlations 
between the two modes. These three real numbers provide an equivalent representation 
of the covariance matrix since inverting Eqs. (jH) and inserting into (jl3p gives 

n + ^ = ^{v^ + r^) , Mc) = l{K-r^), im{c) = r^^. (u) 

As we shall see in Sec. IIIID| the determinant of C is related to S, the von Neumann 
entropy of the state. With (fT2]) . one finds that the x (time) dependence drops out and get 

det(C) = ^ ^ |c|2 = n(n + l) ^ S = 0. (15) 

The minimal value that det(C) can take is 1/4. When this lower bound is saturated, the 
state minimizes the Heisenberg uncertainty relations since det(C) = V^Vt^ — V"^^ > 1/4. 

The state of the linear perturbations is therefore characterized by a single function, 
the power spectrum Vip{k^rj) (and an angle that plays no part in this paper). In realistic 
models, i.e. with interactions, the complete description of the system requires the full hier- 
archy of connected correlations functions. Such knowledge is out-of-reach, and in practice 
one resorts to a coarse grained description. This gives a reduced state prcd characterized 
by a non zero entropy. The choice of the correlations that are discarded must be done on 
physical grounds and in a way consistent with the dynamics. This is the subject of the 
next Section. 

III. COARSE GRAINING AND REDUCED GAUSSIAN DENSITY MATRICES 



In IIII A( we present an operational definition of reduced density matrices appropriate 
for interacting field theories. The advantages of this definition will be emphasized in paper 
II through the analysis of explicit models. We apply it to cosmological perturbations in 
nil CI in the simplest case, the Gaussian approximation. The properties of the Gaussian 
states are summarized in IIII Dl 
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A. Operational definition of coarse graining 



The program begins with the specification of a finite set of "relevant" observables 
(which can be functions of both time and space or momentum). This set 
defines our knowledge about the state and the dynamics (see Appendix |A] where this 
important aspect is emphasized) of the system. In general the finite set 
allows only for a partial reconstruction of this state. The reconstruction is performed in 



the following way. The reduced density matrix pred {^i > 
on the expectation values Oj of these observables, 

Tr (prcd) = 1 , 

TV Led((5i,...,0„)0j 



, On) is defined by the constraints 



(16) 



There is of course an infinite number of density matrices verifying these constrains. How- 
ever, to be consistent with our hypothesis that the Oj are the only observables accessible 
to us, we must choose the density matrix pred which maximizes the entropy given the 
constraints (jl6p : 

(17) 



5" [Pred] > S[p] , 

where 

S[p] = -Tr(plnp) 
The formal solution of this variational problem is [T 



Pred 



exp 




o 



Tr 



exp 



(18) 



(19) 



(In the case these observables depend on space-time, one should read ^Sx X{t, x)C'„(t, x)). 
This is an out-of-equilibrium generalization of Gibb's canonical and grand canonical states 
which reduce to these distributions in equilibrium. 

The Xn are Lagrange multipliers. The constraints on the expectation values ()16p are 
therefore written 



--^In Zo . 



(20) 



One then inverts this system of n equations to express the Lagrange multipliers in terms 
of the expectation values 



A,- 



{Oi,...,On). (21) 

and substitute in (jl9p . The von Neumann entropy of the solution prcd is the Legendre 
transform of the logarithm of the partition function Zp, 



S[p 



red I 



lnZo + ^XjOj = S{Ou...,On) . 



(22) 



For quantum fields, one notices the close resemblance between In Zo and the generating 
functional of Green functions, and between S and the effective action. 

Notice that the solution (llDh is formal. Ambiguities stem from the non-commutativity 
of the operators Oj. As a result, (j22p is strictly valid only when all the Oj commute 
(see also Eq. (IA3P and the following comment). In the following we will only consider 
Gaussian density matrices. The theory of their representation is well developed and we 
shall rely on this body of work. In particular, the ordering ambiguities in (jl9p and (j22p 
can be resolved. The correct formula of the entropy is given by Eqs. (I26p and (ITT 
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B. Application to quantum fields 

For quantum fields, the above program transposes into the following. We start from 
the observation that the knowledge of the Green functions of a (self-) interacting field 
is equivalent to the knowledge of the state of that field. A coarse graining is therefore 
naturally defined by a truncation of the (BBGKY) hierarchy of Green functions at a given 
rank N. As explained above, to be self-consistent, the hierarchy of Green functions must 
be closed. 

This coarse graining is the field theoretic version of Boltzmann's Ansatz. In the latter, 
A^-body interactions are neglected for > 3 and the object of physical interest is the one- 
particle correlation function. Beyond the Gaussian approximation, this coarse graining is 
formalized in terms of the so-called n-Particle Irreducible representations of the effective 



action jl3l ]. 



C. Coarse grained description of metric perturbations 

Let us apply this program to the scalar perturbation ip. The lowest non trivial order of 
truncation is at = 2. The reduced density matrix is then defined by the anticommutator 
of if. The corresponding reduced density matrices pred which maximize S are Gaussian 



14| . These states will be described in details in the next Sections. Here we wish to give 
a qualitative understanding of this coarse graining. 

Each Fourier component of the anticommutator describes the effective evolution of 
the two mode sector ((/9k)<^-k) and can be analyzed separately, as it was the case in the 
linearized treatment. The growth of entropy associated with the coarse graining in the 
present case can be therefore described by the set of Gaussian two-mode density matrices 
/^iT-k- Each of them characterizes the loss of entanglement between the two modes in the 
presence of interactions, which are self-interactions of gravity or (and) interactions of the 
metric perturbations with other fields (e.g. the fields of the Standard Model of particle 
physics or isocurvature perturbations in multi-inflaton field scenarios) . The environmental 
degrees of freedom are thus either the collection of the modes (pk'^^k, or modes of other 
fields. 

Returning to the formalism, for each two-mode, pj^^-k defined by the anticommutator 
Gaiv,v';k) = \ ld^xe'''^TT[ptot{mn) {v^(r?,x),v9(r/',0)}] , (23) 

where the in-state ptotiVin) is the interacting vacuum of the total system. It replaces 
the Bunch-Davis vacuum of linearized modes. ^p{r], x) is the Heisenberg operator of the 
nonlinear metric perturbations. In the SP, using (jlip . Prcd(??) is the Gaussian density 
matrix possessing the following covariance matrix 

Vip = Gair], rj; k) , V^^ = dr,'Gair], rj'; k)[rj=ri' , 'Pn = dndr,'Ga{r], f]' \ k)U=-n' ■ (24) 

It will indeed be shown in [3] that in the Gaussian approximation it is always possible 
to make a canonical transformation {f,TT') (<^, tt = dy(p). As recalled in Sec. IIVI the 
entropy is invariant under canonical transformations. Throughout this paper we use this 
Gaussian approximation. Since the interactions do not spoil the property of statistical 
homogeneity, the reduced states can be described in the same way as the pure state of 
the linear perturbations [Sec. HI]. Only the values of the covariance matrix element differ. 
These are now given by (j24p in place of (jl2p . We call these states Gaussian Homogeneous 
Density Matrices (GHDM). 
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D. Entropy of picd 

To measure the strength of the correlations between k and — k, we introduce the pa- 
rameter 6 defined by 

|cp = + 1 - 5) , 0<(5<n+l. (25) 

The standard inflationary distribution of Section II is maximally coherent and corresponds 
to 5 = 0. The least coherent distribution corresponds to 5 = n + 1. The entropy (jl7p is a 
strictly growing function of 6. It is related to the determinant of the covariance matrix by 

S = -Tr (plnp) = 2 [(n + 1) ln(n + 1) - nln{n)] , (26) 

where the parameter n is defined by 

. 1 



n+-j =det{C) = - + n5. (27) 

The prefactor 2 in (|26p accounts for the fact that p is the state of two modes. In the range 
6 ^ 1/n, n ^ 1, the expression (|26p of the entropy simplifies, 

n~V^, S = ln{n6) + 0{l). (28) 

These equations have a simple geometric interpretation that we will use in Sec. IVII and 
IVII Al The instantaneous eigenvectors of the covariance matrix are given by the rotated 
canonical variables 

$k(7?) = COS e Vk^^ + sme ^ = — (ajf e"^^ + a*_"Je+'^) , 

nk(7?) = - sin Vk^p'^^ + cos 6* ^ = ^{a^e-'^ - a^^^e+'^) . (29) 

The angle 6 = 9{r]) = ^ arg(c(?7)) gives the orientation of the eigenbasis of the covariance 
matrix with respect to the (fixed) original variables ip^ and vr^^j^. The eigenvalues of the 
covariance matrix are the variances of $ and 11, 

({$k, 4}) =n + ^ + \c\ = 2n + 0(1, 6) , (30a) 
({Hk, 4}) = n + i - |c| = ^ + O (i) , (30b) 

and ({'I'k) n|^}) = 0. Hence, when 6 ^ 1/n, one has 

det(C) ^nS^A =^ S^lnA, (31) 
where A is the area under the 1—a contour in phase space. 



IV. THE INTRINSIC PROPERTIES OF GHDM 

The von Neumann entropy S = — Tr (pln(p)) is manifestely an intrinsic property of the 
state. We show that for GHMD it is the only intrinsic property. Important consequences 
are then derived. 
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A. Entropy as the unique intrinsic property 

To illustrate the question, let us first consider the following situation. In place of the 
Mukhanov- Sasaki variable, one could choose instead to work directly with the curvature 
perturbation C. The quadratic part of the Lagrangian is 



c 



^ jdtd'xe{t)[a\C?-a{VC?], (32) 



where C, = dtC = o ^^nC- The conjugate momentum of the curvature perturbation is 



^cW = 4;^C(t). (33) 

Assuming the slow roll condition e ~ cte and e/He ~ cte, the modes with positive fre- 
quency in the asymptotic past are 

Ckiv) = C (1 + ikr^) e-''' , ICT = ^ ^ • (34) 

From this solution and (|33|) , one obtains the following expressions for the covariance matrix 
at tree level 

n = ^({Ck(r),C-k(r)})=|C,°|'(l + x2) , (35a) 

= ^({Ck(r), vr_k(r)}) = Hx'\C\' , (35b) 

^ ^({vrk(r), vr_k(r)}) = H^x^'.f . (35c) 

Despite the differences with the variances (|12p . one checks that det(C) = 1/4. The 
determinant of the covariance matrix, hence the entropy, are invariant under the linear 
canonical transformation (y?, vr) (Cvr^^). 

More generally, the entropy does not depend on the choice of canonical variables. 
Indeed, canonical transformations are linear symplectic transformations of the covariance 
matrix. The intrinsic properties of /Ored are therefore the symplectic invariants (of the 
group Sp(4, R) for a system of two modes). The symplectic invariants of a 4 x 4 covariance 
matrices are known [l5|. When the constraint of statistical homogeneity is added, these 
invariants are degenerate into a single quantity, namely the determinant of C. In other 
words, det(C) is the unique intrinsic property of the effective Gaussian state pred of the 
cosmological perturbations. Unique here is employed in the sense of an equivalence class: 
any quantity which can be expressed in terms of det(C) only, e.g. the entropy, characterizes 
the same quantity. It is interesting to notice that this uniqueness rests on statistical 
homogeneity. 

We now apply this result to prove that the entropy is well defined despite the redun- 
dancy of Einstein's equations and the arbitrariness in the definition of the perturbatively 
renormalized Green functions. 



B. Entropy and gauge invariance 



The entropy of metric perturbations is independent of the choice of gauge. The reason 
is that in a change of coordinates + S,^, the Lagrangian (of Einstein-Hilbert plus 
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inflaton) changes by a total derivative. It can therefore be seen as a canonical transforma- 
tion. That the entropy is a gauge invariant quantity can also be seen independently from 
the fact that we quantize gauge invariant variables like the Mukhanov-Sasaki variable or 
the comoving curvature perturbations. 



C. Entropy and renormalization 

To show the independence of the entropy of the renormalization scheme, we proceed 
in three steps. We consider first renormalizable field theories in Minkowski space, then 
generalize to non renormalizable theories in Minkowski space, and finally to non renor- 
malizable theories on a curved background. In what follows we are only considering the 
entropy per two-mode. 

For renormalizable theories in Minkowski space, we have to examine two potential 
sources of ambiguities, the parameters of the Lagrangian and the wave function renormal- 
ization. First consider the parameters. They are local terms in the effective action, and 
therefore appear as local terms in the equation of the propagator. It can be checked (see 
paper II) that they are not responsible for the change of entropy. 

In a renormalizable theory, after renormalization of the bare parameters, the renor- 
malized Green's functions are numerically equal to the bare ones, up to a multiplication 
factor given by the wave function renormalization constant(s) to the appropriate power. 
The field strength renormalization is also a canonical transformation, and therefore leaves 
the entropy unchanged. The quadratic part of the bare Hamiltonian is in general 

1 r / 7r2 



H = -jd'x [Z.^^ + Z(^K'ea{VCYj , (36) 

where Zj^ and renormalize the operators (7r(x)7r(?/)) and {(^{x)(^{y)) . The commutation 
relations [C(i,x), 7r(t,y)] = i6^^^ (x — y) implies 

Zn = ^. (37) 

The renormalization of the wave function therefore leaves det(C) unchanged since the two 
products (Ck(i)C-k(i)) x (7rk(i)7r-k(i)) and {{Ck{t), 7r_k(i)})^ involve the same number of 
C and vr. 

The situation is a little different in a non renormalizable theory (still in Minkowski 
space) because the cancellation of divergences at the order L of the loop expansion requires 
counterterms with higher derivatives. However at order L, these counterterms are local 
and therefore do not change the entropy. In other words, the non-local part of the effective 
action (e.g. in Fourier space terms like ln(— fc^/ /i^) or y^^^W/Jj?) responsible for the growth 



of entropy are precisely those predicted by the quantum theory [16|]. Irrelevant operators 
containing more that two derivative require some care, since the equation of the two-point 
function is then higher than second order and can no longer be solved from the knowledge 
of the covariance matrix at an initial time. These higher derivative terms are generally 
discarded in a self-consistent perturbative treatment, because the extra solutions of the 
propagation equation are not analytic in the limit h ^ In this self-consistent 

perturbative sense, these counterterms do not introduce ambiguities to the definition of 
the entropy. 

Transposition to a curved background should not alter these conclusions because the 
previous considerations about the canonical structure, the distinction between local and 
non-local corrections and their distinct contributions to the entropy, and the role of higher 
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derivative operators, still hold on a curved background. However, we have not verified 
explicitly for the models of paper II that the counterterms do not introduce ambiguities 
in a self-consistent perturbative treatment. We leave this analysis for future work. 



V. SEPARABLE DENSITY MATRICES AS CLASSICAL STATES 

To interprete the value of the entropy induced by a given coarse graining, we need 
a classification of the corresponding reduced states. In the class of two-mode density 
matrices there exists an operational definition of classical states. After recalling the defi- 
nition of separability, we explain why it cannot be applied unambiguously to cosmological 
perturbations during the inflationary era. 



A. Definition and characterization of separable states 



The state @ is entangled, i.e. it violates some Bell inequalities [18|]. We recall that 
this entanglement refers to the correlations between the modes of opposite wave- vectors. 
The definition of classicality called separability is a statement on the nature of these 
correlations. Explicitely, separable states do not violate any Bell inequality based on pairs 
of observables, each of them acting only in one sector, i.e. on k or on — k. These states 



can therefore be represented as a convex sum of tensor products of density matrices [19|], 

Psep = ^Pipi^ ^P%, Pl>0. (38) 

I 

Unlike entangled states, these can be prepared with local classical operations (in the 
sense that one can prepare the states p^ and p^^^ without making the modes k and 
— k interact), and by using a random generator (characterized by the probabilities pi). 
Although physically transparent, this definition (j38p is of little practical use because of 
the difficulty of proving or disproving the existence of the set of p^^^ and pi. Fortunately, for 
Gaussian states, a criterion of separability is known [20I ]. For the GHDM, it is expressed 
as the following inequality on the parameter S 0] 

Pred separable <^=^> 6 = 1. (39) 

In brief, GHDM fall into two disjoint classes, separable states {6 > 1) which are opera- 
tionally indistinguishable from stochastic ensembles, and entangled states {6 < 1) char- 
acterized by a departure of the anticommutator from classical correlations. This is the 
definition of a classical GHDM we adopt. 

Finally, at the threshold of separability, the entropy is equal to one half of the entropy 
of the thermal state with the same value of n (for n ^ 1), see Eq. 



5(n, 6 = 1) = ^SmM ^ ln(n) . (40) 



B. Discussion 



From (^^, we see that the values of n and c, and therefore that of 6 as well, depend of 
the choice of canonical variables. As explained in Section [TVl only the combination 

det{C) = (n + ^y -\c\\ (41) 
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does not. This poses a fundamental limitation to the applicability of the criterion of sepa- 
rability to systems with a time dependent Hamiltonian, and therefore to the cosmological 
pertm'bations during inflation, see ([2|). 

For instance, in place of the variables (9?k)7'"-k = dr^fk) and the Hamiltonian ([2]), one 
could choose to work with the variables (p]^ and -7f_k = Sj^V'k ~ {z' /z)(pi^, as done for 
instance in Q]- In this case the corresponding Hamiltonian is 

d z 

^k,-k = Ivi-kP + ^^Iv^kP + — (TT-kV'k + TTkV'-k) • (42) 

z 

We leave it to the reader to verify that in the Bunch-Davis vacuum the functional de- 
pendence of the variances is different from Eq. (jl4p while det(C) = 1/4. However, the 
canonical transformation from these two sets of variables mixes modes of opposite wave 
vectors, i.e. it belongs to Sp(4) but not to its subgroup Sp(2) x Sp(2). The property of 
(k, — k)-separability is therefore not stable under such canonical transformations. 

As a result, the threshold of separability (j40p depends, during inflation, on the choice of 
canonical variables. The inequality 5 > 1 is in fact an inequality on the eigenvalues of the 



so-called partial transpose of the density matrix [20(] . Namely, the lowest eigenvalue should 
be larger than the variance in the vacuum state. Hence, for any choice of the creation 
and annihilation operators (a(j),o|^p, it yields > 1. Therefore a separable state in 
i-representation can be entangled in other representations related to it by a Bogoliubov 
transformation mixing the modes k and — k. 

In this we find another illustration of the well known fact that when the frequency 
varies, one looses some of the usefull characterizations of quantum states. Remember that 
in time dependent backgrounds there is no intrinsic definition of the occupation number, 
even though the expectation value of the stress tensor stays well defined. Here, even 
though the entropy is well defined, the criterion aofclassicality is not clearly defined as 
long as the frequency significantly varies (to be more precise, outside the validity domain 
of the WKB approximation.) Since there is no clear notion of particle as pair creation (or 
mode amplification) proceeds, there is no clear distinction between quantum and classical 
states, as illustrated by the criterion of separability. 

On the contrary, during the radiation dominated era, this ambiguity disappears since 



the mode frequency of the Mukhanov- Sasaki variable is constant (because z = cte) |lll ]. 
Then, defining classical states as the separable ones presents several advantages. First, 
since it relies on the possibility of violating Bell inequalities, it is an operational definition. 
Second, the separation between quantum and classical states is sharp. This contradicts the 
common belief that the quantum-to-classical transition is fuzzy. Being sharp, "the time 
of decoherence" is also precisely defined. Starting from an entangled state, this transition 
occurs at the time when 5 crosses 1. 



VI. OTHER CRITERIA OF CLASSICALITY 



We show the equivalence between separability and three other classicality criteria. The 
latter suffer from the same ambiguities during inflation as the one of separability. They 
must therefore be compared to each other in the same representation, e.g. that defined 
in Sec. [Hi We conclude by discussing a class of inequivalent criteria which do not suffer 
from the above ambiguity but which are based on another arbitrary choice. 
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A. Criteria equivalent to separability 

1. The broadness oj the Wigner representation 



A first alternative criterion was introduced in [21[ for one-mode systems, generalized in 
[i^ l to two-mode systems, and applied to cosmological perturbations in IC]. It rests on the 
observation that in many respects, the Wigner representation behaves like a probability 
density over phase space, except for the fact that it can have a finer structure than h, and 
in particular can take negative values in small regions. 

We recall that the Wigner function can be defined by (for a system with canonical 
variables {q,p)) 

e-^%(,-f,,4), (43) 

For Gaussian states, Wp is Gaussian and its covariance matrix is C of Eq. (jlip . Although 
it is positive everywhere, we see from Eq. (|30b ) that when 6 < 1, the variance of the 
subfluctuant variable 11 is smaller than the variance in the vacuum. The criterion of clas- 
sicality therefore consists to ask that the Wigner representation does not contain features 
that are smaller than those of the Wigner function in the vacuum state (this definition 
includes the case of negative values). For GHDM, we therefore expect that this is the case 
when 

(nnt) > 1 , (44) 

i.e. 5 > 1 + 0(1 /n), which is equivalent to the criterion (|39|) for large n. 

This requirement is made more precise by asking that the Wigner representation of prcd 
is broad enough to be the Husimi (or Q-) representation of some normalizable density ma- 
trix p'. Indeed, the Wigner representation of a state p' is mapped onto its Q-representation 
by a convolution with a Gaussian function of covariance matrix ^, i.e. 

QA<1,P) = /^e-(-^')^-(--')V,,(<z',p') , (45) 

which explains why the Husimi-representation Qp' of p' is a broader function than its 
Wigner representation Wp/. Moreover, the Husimi-representation Qp/ of any density 
matrix is positive because Qp'{q,p) is the expectation value of p' in the coherent state 
\{q + ip)/\/^). We show in Appendix |B] that the GHDM p verifying this condition are 



5>5Q = l + i^=4ep + ^, (46) 



as anticipated from the heuristic argument in the previous paragraph. It is larger than 
the condition of separability, but only slightly, and the two criteria are equivalent in the 
limit n ^ 1 relevant for cosmological perturbations. 



2. P-representability 

A second alternative definition of classicality is the requirement that pj-d admits a P- 
representation as defined by Glauber [1, [2^. It means that the states can be written as a 
statistical mixture of coherent states, 

P{v,w)\v,'k){v,'k\ (g) \w,-k){w,-'k\ . (47) 

vr vr 
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where P{v, w) is a normalizable Gaussian distribution. A P-representable state is obvi- 
ously separable. In general, the converse is not true but it turns out that for GHDM 
P-representability and separability are equivalent. Hence a GHDM is P-representable if, 
and only if 

5 > 1 = 5scp . (48) 

To avoid any misunderstanding, we emphasize that this condition does not mean that 
coherent states are "the pointer states". These have only been used as a resolution of the 
identity. We shall return to this point in the next Section. 



3. Decay of off- diagonal matrix elements 

A third alternative condition of classical two-mode states is provided by the decay of 
interference terms of macroscopically distinct states, sometimes referred as Schrodinger 
cat states. Since this decay is asymptotic, this criterion is less precise then the other three 
criteria. However it is qualitatively equivalent to them in the following sense. We refer 
to Appendix D in 0] for details. The off-diagonal matrix elements of the density matrix 
of the pure state ([6]) in, say the basis of coherent states are correlated over a range oc n. 
Using the coherent states as representative of semi-classical configurations of the field 
at a given time, the squeezed state is a linear superposition of macroscopically distinct 
semi-classical configurations. As for the entropy, see ()40p . the correlation length between 
off-diagonal matrix elements is very sensitive to the value of 5 in the range [0, 1] where 
it decreases monotonously from 0(n) to 0(1) as 6 increases. For 5 > 1, the correlation 
length depends very slowly on 6 and stays 0(1). Hence the decay of the correlations length 
also distinguishes classical Gaussian states as those with 6 > 1. 



4-. Adding one quantum incoherently 
There is a simple physical interpretation to the criterion 6 > 1. Such a density matrix 



is obtained from the pure state by adding incoherently one quantum on average 23| . One 
obtains 

n t-^ n' = n + - , |c|i— >|c'| = |c| = n{n + 1) . (49) 

We have split evenly the contribution of the quantum between each mode in order to 
preserve statistical homogeneity. One obtains that the corresponding value of 5 defined 
by |c'p = n'(n' + 1 - 5) is 

5 = l + ^ + o(^, (50) 

in agreement with the other criteria when n ^> 1. As an interesting side remark, re- 
membering Eq. (I40p . the fact that the entropy gain associated with this addition of one 
quantum is large clearly establishes the fragile character of the quantum entanglement. 



B. Inequivalent criteria 



Let us now discuss a criterion [10|, |2j| which is not equivalent to the above four criteria. 
Instead of adding incoherently one quantum, one can consider the entropy gain associated 
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with the loss of one bit of information. Then, whatever the system is, the change of entropy 
is 5 = ln(2). Applying this criterion to the standard inflationary distribution, the value 
of the parameter 5 is 

(Joncbit - — ■ (51) 
n 

The statistical properties of the corresponding density matrix are essentially the same as 
those of the original pure state ^. For instance. Bell inequalities are still violated and the 
correlation length between off-diagonal matrix elements is still 0{n). 



C. Summary 



In conclusion, we have found that the various criteria of decoherence fall into two 
distinct classes. On the one hand, classical states with respect to the statistics of the 
anticommutator O = {(pk{t), separable states with 6 > 1. This criterion 

is ambiguous during inflation because it rests on a choice of canonical variables ^p^jTr^. 
On the other hand, a lower bound on the entropy S > ln(A'^) is intrinsic as only the value 
of the entropy is involved. However, we could not identify any operator(s) to which this 
criterion might refer to. The value of is therefore not dictated by any physical property 
of the state. (It might be provided by the resolution of observational data, but this would 
confirm that it does not characterize the state of the system.) 



VII. DECOHERENCE AND STATISTICAL MIXTURES 



The two classes of criteria presented above seem a priori rather different as the first 
class is based on the quantum properties of the system, whereas the second class uses 
only the value of the entropy. In this Section, we show that these two criteria can be 
incorporated into a single treatment. 



A. A picture of decoherence 



We need first to express two-mode density matrices has a tensor product of two one- 
mode density matrices. Indeed, recall that entangled states cannot be represented in the 
form (j47|) . Instead all GHMD can be decomposed as 

Pk,--k{S) = pi{S)®P2{S), (52) 

where 1 and 2 refer to a separation of the Hilbert space into two sectors defined by the 
variables 

V^i,2 = • (53) 



^ It was noticed in [U that adding one quantum yields to S = S'max/2, but its physical interpretation 
was not mentionned. In this early work, the authors do not however refer to criteria of section IVI Al 
and prefer instead the criterion S = ln(2), referring to [i^ for an experimental situation where deco- 
herence is effective for such values. This reference is misleading. Indeed, the experimental observation 
of decoherence reported there concerns a system of a two-level atom in a cavity interacting with the 
electromagnetic field in a coherent state (the environment). One finds that the interferences are blurred 
when one photon is exchanged. Since the Hilbert space of the atom has two dimensions, the exchange 
of one photon between the system and the environment is equivalent to the exchange of one bit. In 
inflationary cosmology since n ~ lO'^"", this correspondence is lost. 
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Because of homogeneity, the matrices pi and p2 are characterized by the same covariance 
matrix C, which moreover coincides with that of p. Exphcitely, one has 

C, = C,= ( ] = (l- l--)=C. (54) 

In this way, the properties of the state of cosmological perturbations have been encoded 
into two fictitious one-mode systems. The entanglement {6 < 1) between modes of opposite 
wave vectors k and — k reflects into the existence of two sub-fluctuant variables Ili^2 as in 
Eq. ([sob). 

The question we address concerns the use of minimal Gaussian states \{v,^)) of V'i,2) 
i.e. squeezed coherent states, to represent the states /9i,2(^) as statistical mixtures in the 
following sense 

Pi(^) = I^Pd^;6)\{v,0){{v,0\, (55) 

that is, by summing only over the complex displacement v. The basis of states we use 
therefore have a common orientation and elongation which is fixed by the squeezing pa- 
rameter 

^ = re^'^^. (56) 

Using the la contour in phase space, the state \{v,(,)) draws an ellipse of unit area, 
centered around {(pi cx Re(z;), vfi oc Im(t;)), with a long axis (the superfluctuant direction) 
oc e^^ making and angle 6c w.r.t. the horizontal ipi-axis. The latter is chosen for 
simplicity along the superfluctuant mode of pi deflned at Eq. (j30b ). Hence in this "frame" , 
9c is the relative angle between the big axis of pi and the big axis of |(w, ^)), see Eq. (IC12p . 

As shown in Appendix[Cl the states pi,2{S) can be represented as in ([55]) for any value 
of 6. More interestingly, the choice of the basis vectors, which is parameterized by ^, is 
more limited when (5 < 1 in that the range of allowed values of Oc belongs to a bounded 
interval which shrinks to zero as 6 ^ 0. That is to say, the pure state {5 = 0) admits 
only one representation, itself. As 5 increases from zero, the allowed range of 9c increases 
but is necessarily strictly smaller than 7r/2. One can say that pi "polarizes" the pavement 
()55p along its superfluctuant mode. In addition, the range of r also increases and the 
distribution P^{v) becomes broader, which means that a growing number of families of 
minimal states can be used to represent pi{S). For any 6 < 1, one must have r > 0, i.e. 
the states \{v,S,)) are necessarily squeezed (and as we saw they tend to align with the big 
axis of pi). When the threshold of separability is approached, i.e. 5^1", the lower 
bound of r decreases to zero. In this limit one can represent pi{5 > 1) with coherent 
states. Notice also that as 5 increases, the distribution V^{v;6) becomes broader (in v) 
and at the threshold of separability does not have any structure finer that a unit cell of 
phase-space. In this we recover what was observed in subsection IVI A Tl 

The fact that pi polarizes phase space is easy to understand a contrario. Indeed, 
consider the minimal states which are squeezed in the direction perpendicular to that of 
pi, i.e. 9c = 7r/2. The corresponding spread in the subfiuctuant variable 11 of (129]) is 
large. A statistical mixture of these states is necessarily spreaded out in the 11 direction, 
i.e. (nnt) > 1/2, that is (5 > 1, see dSO]), and therefore S > S'Max/2. 

Let us formulate these results the other way around in terms of 6^ (or S^), the amount 
of decoherence (or entropy) needed for the state to be written in the ^-basis as in Eq. (|55p . 
According to what we just said, the further |(v,^)) departs from pi, the larger is 6^ (and 
5g). For well aligned |(f ,^)), i.e. for 0c = 0, V |^| = r ^ 0, we have 5^ < 5scpi the minimal 
amount of decoherence to reach the separability threshold. 
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Hence given a certain decoherence rate, these considerations translate into the time 
of decoherence t^, i.e. the time after which any initial state /Ji(to) has evolved, as 
decoherence proceeds (as 6 increases), into the statistical mixture ([SSj) for that value of ^. 
When using again well aligned \{v,$,)), this defines ts^p as the maximal t^. It is also also 
the first time such that the representation (j47p is allowed. 



B. Ambiguity in choosing pointer states 

Pointer states are meant to bridge the gap between the quantum and classical descrip- 



tions of a system [271 . 1281 ] . The defining property of pointer states is their robustness over 
a given lapse of time tp.s. which makes them the quantum counterparts of points in the 
phase space of classical mechanics. That is, once the system has been prepared into a 
given initial state and placed into contact with a given environment, they are the states 
the least perturbed by the environement over tp,s.- One generally obtains radically dif- 
ferent pointer states whether tp,s. is much smaller or comparable to the dynamical time 
scales of the open dynamics (the proper frequency of the system or the characteristic time 
of dissipation). 

To make the choice of ip.s. less arbitrary, one often adds the requirement that any 
initial state of the system evolves, over a time to, into a density matrix which cannot 
be operationally distinguished from a statistical mixture of the pointer states. Since the 
evolution of both the pointer states and p are governed by the same dynamics, consistency 
requires that the times tp.s. and tu he commensurable. 

As we saw at the end of Sec. IVII A\ choosing a family of pointer states (i.e. a pair 
{r,6c)), is arbitrary since it amounts to a choice of 5^, or to a choice of a lower bound 
S > ln(A''g) which we recall is not dictated by physical considerations. The corresponding 
time of decoherence t£) = is therefore arbitrary. But there is a more fundamental 
obstruction to the identification of a pointer basis during inflation, namely pointer states 
refer to a choice of canonical variables, which is arbitrary when the Hamiltonian depends 
explicitly on time (compare for instance narrow wave-packets in or C). 

In brief, on the one hand during inflation, the question of finding "the" pointer states 
is not well-defined. On the other hand, during the radiation dominated era, the criterium 
of separability offers an unambiguous definition of classical states based on the statistical 
properties of the state. From these two facts we conclude that the concept of pointer 
states does not seem useful to analyse the decoherence of cosmological perturbations. 



VIII. SUMMARY 

By truncating the hierarchy of Green functions, we first show how to get a reduced 
density matrix for the adiabatic perturbations, in a self-consistent manner, and from the 
interacting properties of the system itself, i.e. without introducing some ad hoc environ- 
mental degrees of freedom. 

When truncating the hierarchy at the first non trivial level, statistical homogeneity still 
implies that the density matrix factorizes into sectors of opposite wave vectors. Hence the 
reduced density matrix of each sector is determined by three moments related to the anti- 
commutator function, see ()24p . This also implies that decoherence here describes the loss 
of the entanglement of these two modes. The level of decoherence is characterized by one 
real parameter in each sector, the parameter 5 introduced in ()25|) . 

We then show that the entropy S contained in each reduced two-mode density matrix 
is a well-defined quantity which monotonously grows with 6. The important conclusion is 
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that S is the only intrinsic quantity of these reduced density matrices, in that the other 
quantities ah require to have chosen some pair of canonical variables to be evaluated. 

After the entropy, we studied the quantum-to-classical transition. We show that the 
criterion of separability agrees with three other criteria, namely the broadness of the 
Wigner function, the P-representability, and the neglect of off-diagonal elements of the 
density matrix. During inflation, these four concepts are ill defined, because, contrary to 
the entropy, they require to have selected some creation and destruction operators, which 
is ambiguous since the frequency significantly varies when modes are amplified. 

We compare the above four criteria which give a large entropy at the threshold of 
classicality, see ()40p . to another class of criteria which give much smaller entropies, see 
(j5ip . and explain why they differ so much. In the last Section, we present a unified 
treatment of the reduced density matrices in which both types of criteria can be found. 
We show that each reduced density matrix can be written as statistical mixtures of minimal 
states which posses a well defined range of the squeezing parameter ^ of Eq. (f56l) . The 
details are given in Appendix [Cl 

This analysis clearly shows that there is no intrinsic definition of a threshold of deco- 
herence (or a critical entropy) at which the quantum-to-classical transition would occur 
during inflation. We also argue that the pointer states are of no help in providing such a 
definition. 

In the next paper, we calculate the entropy for two dynamical models. It appears 
that no significant entropy is gained during single field inflation, so that the quantum- 
to-classical transition should occur during the adiabatic era. On the contrary, it is very 
efficient in multifield scenarios. 



The work of D.C. is supported by the Alfried Krupp Prize for Young University Teachers 
of the Alfried Krupp von Bohlen und Halbach Foundation. 

APPENDIX A: WHY BOTHER ABOUT THE COVARIANCE MATRIX ? 

In this Appendix, we emphasize the role of the dynamics and the canonical structure 
in defining reduced states and their associated entropy. 

One may object to our analysis of the two-point function at different times on the 
grounds that one can only measure the values of this correlation function on our past 
light cone or on the last scattering surface. This objection can be answered as follows. 
The outcomes of a measurement are analyzed through the grid of a particular dynamical 
model. There are in general more dynamical variables necessary to write a consistent 
model than one can actually measure. The state of the system, either the exact state 
or a partial reconstruction of it, depends on both the outcome of the measurement (for 
instance, the values of the power spectrum and bispectrum) and the correlation functions. 

As an example, we keep only the power spectrum and show that this assumption is 
untenable as it amounts to ignore the canonical structure of the theory. As a result the 
von Neumann entropy of the reconstructed state is not well defined. Indeed, following the 
algorithm of Sec. IIII Al we write the Ansatz for the density matrix 



Acknowledgments 




) 





(Al) 



19 



where A is a Lagrange multiplier ensuring that the power spectrum of C has the measured 
value P^, i.e. V^iq) = Tr (/OcCqC-q) = 1/2A. The distribution P{C) is therefore 



and the von Neumann entropy is then found to be 



S = \nZ^ + XV^ = ^ In (27rP^) + ^ . (A3) 

The constant 1/2 is universal and corresponds in this scheme to the entropy of the vacuum. 
For ~ lO"^'^, the entropy (|A3p is negative. To cure these pathologies, we must add 
some physical input, i.e. enlarge the set of observables to include the variances Vn and 

Another way to see the necessity to consider C rather than Vf^ alone is the following. 
Note that even though the state ()A1|) is Gaussian, we cannot calculate its von Neumann 
entropy from the formulas (jlip . (j26p and (|27p . The latter are only valid for a state 
reconstructed from the anticommutator function. If one insists on doing so, one finds 
Vtt = oo (since is diagonal in the field-amplitude basis), so that the entropy is infinite. 
The state must be therefore be regularized first, by giving a finite width to (vr^). This 
operation is arbitrary without additional physical input about the dynamics. 



APPENDIX B: EQUIVALENCE OF THE CRITERIA OF SEPARABILITY AND 
OF BROADNESS OF THE WIGNER FUNCTION 



Let us consider a Gaussian density matrix p of a bipartite system and let C be its 
covariance matrix. We adopt a different parameterization as in the text, following 20|] 



C = Tr(p [a, At}) 
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The positivity of the density matrix and the non-commutativity of the creation and anni- 
hilation operators puts a constraint on C, namely 



E 

C+->0, 
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This is the general form of the Heisenberg uncertainty relation for the second moments. 
Indeed, this condition puts a lower bound on the lowest eigenvalue of C, which, in the 

case of homogeneous states, reads n + 1/2 — \J\c\^ + 1/4 > 0. The latter can be recast as 
the Heisenberg uncertainty relation 



(aka[.)(aj^ak) > |(akak>|' 



(B3) 
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A necessary and sufficient condition for the separability of a GHDM ^ is 



E 

p separable <^=^ AC A + — > , A 
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(B4) 



which gives n — |c| > 0, or (5 > 1 using (j25p . 

A different criterion than separability is used in 211]. With this condition, the Wigner 
representation of p must be broad enough to be also the Q-representation of a density 
matrix p' . We specialize to Gaussian states. To obtain the corresponding condition on 
the Wigner function, one simply notices that the convolution of the Wigner representation 
of a state p with a Gaussian of covariance matrix 1/2 yields the Q-representation of that 
state. 
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where = (93, vr^, 99^, tt) . Hence the condition ()B2|) on C can be written as a similar 
condition for the covariances Cq of the Q-representation 



C + \{E-l)>Q. 



(B6) 



In consequence, a necessary condition for the Wigner representation Wp of a state p to 
be also the Q-representation Qpi of a state p' is that the covariance matrix C of p verifies 
C+^(£' — 1) >0. For Gaussian states, this condition is also sufficient. Specializing now 
to the case of GHDM, we arrive at the conclusion 
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(B7) 



This is slightly more constraining than the separability condition n > |c|, but the difference 
between the two is not relevant when n ^ 1. 



APPENDIX C: REPRESENTATIONS OF PARTIALLY DECOHERED 

DISTRIBUTIONS 

In this Appendix we show that all partially decohered Gaussian distributions can be 
written as statistical mixtures of minimal states which belong to a certain family. As 
decoherence increases, the ranges of the parameters characterizing this family become 
larger. As the threshold of separability, the angle between the super-fluctuant modes of 
the minimal states and that of the distribution becomes unrestricted. 



We recall that a necessary condition for p to be a separable density matrix is that its partial trans- 
pose is a bona fide density matrix [2^. The partial transpose ppt is by definition, obtained from 
p by a transposition in one sector only, say — k. In any basis jn,k)|m, — k), this is written as 
(n,k|(m, — k|ppt|m', — k)|n',k) = {n,k|(m', — k|p|m, — k)|n', k). In Eq. (|B4|) this operation is im- 
plemented by the matrix A. Gaussian states are separable when this condition is satisfied [20| . 
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The statistical homogeneity allows a formal reduction of the problem. One can decom- 
pose the field amplitude 4>q into its "real" and "imaginary" parts (f)q = ((/>i + i(f>2)/V^, 
such that 

1 —i 
ai = -j= (ak + a_k) , 02 = (ok - O-k) ■ (CI) 

With this decomposition of the Hilbert space, GHD factorize 

Pk,--k = Pi^ P2, (C2) 
where pi = p2- In addition, the parameters of (jl3p are given by 

n = Tr ^pi a\ai^ , c = Tr (pi a^) . (C3) 

Similarly, for two-mode coherent states \v,]i)®\w, —k) = \vi)0\v2) where v = {vi+iv2)/V^ 
and w = {vi — iv2)/V^- 

Let us consider a Gaussian density matrix p of a single mode, characterized by a value 
ofO<(5<n + l. We ask whether there exists a family of minimal Gaussian states | (f , 0) 
and a Gaussian distribution P(^{v) such that 

PS = I^Pdv)\{v,0){{v,0\- (C4) 
The minimal Gaussian states are chosen to be the displaced squeezed states 

\{v,0)=D{v)S{0\0), (C5) 
where D the displacement operator and 

5(e) = ei(^""''-^-""'), (C6) 



the squeezing operator. The distribution P^{v) is centered and is therefore defined by its 
covariance matrix 



P,iv) = (detQ)-/2 exp {-IxtC-x} , ^ = ( ) , = ( ? n J 



(C7) 



The moments of the state p defined at Eq. (I13p are obtained from the ones of the distri- 
bution Pf by the definition (jCSP 



c = I^Pdv){{v,0\a'\iv,0) = j'^PM + ■ (C8) 



Here, a and (3 are the Bogoliubov coefficients associated with the squeezed state 1^) by 

a = ch(r), /? = e-^^2^=sh(r) , ^ = ^e^^' • (C9) 
Hence, the momenta of the distribution P^ are 

n = n5 + |/3|2, c = Ci: + af3. (CIO) 



22 



The only constraints on and arise from the fact that the right hand side of ()C4p 
must be positive and normahzable in order to be a density matrix. These conditions are 
respectively 

p > -P? > , 

Tr(p) = l ^ J^p^(v) = l ^ det(Q)>0 (Cll) 

We write 

where 6c is the angle made by the squeezed coherent states with the super-fluctuant mode 
(the eigenvector of C with the largest eigenvalue) . We now look for the range of values of 
X = e^^ and 9c allowed by the constraint on the determinant (jClip . that is 

R{x)=Ax'^ + 2BX + C <0, (C13a) 
A = n + ^-\c\cos{2ec) , B = -(^n6 + ^^, C = n + ^ + \c\ cos{29c) . (C13b) 

Since A > 0, the inequality (lC13b ) can only be satisfied if the discriminant C = — AC is 
positive. Since the coefficients A, B and C depend only on 9c and 5, this gives an implicit 
constrain equation for 9c{Si), 

C{9c, 6) = {n6f - \cf sm^{29c) > ^ sm\29c) < g{6) ^ . ^""^ ^. . (C14) 

n\n + 1 — oj 

The function g(S) is strictly growing over the interval < (5 < n + 1 and takes the special 
values (jr(0) = and g(\) = 1- We distinguish the two following cases: 

1) If (5 > 1 (the two-mode state from which p is obtained is separable), all the values of 9c 
are allowed and the squeezed coherent states in ()C4p can have an arbitrary orientation. 

2) If < (5 < 1, the angle 9 can only vary into the interval [—9m, +^a/] where the angular 
opening is defined by 

sm^(2gM) ^ , ■ (C15) 

n[n + 1 — 0) 

The value 9m = 7r/2 corresponds to (5 = 1. The interval [—9m ■,+9m] shrinks to zero as 5 
decreases. In the limit 5 = (that is, for the pure squeezed state), the squeezed coherent 
states 1(^,0) ™ust be aligned with the superfluctuant mode of pi. 

We now specialize to the case < 5 < 1 and we look for the range of allowed values of 
X for a given 6 and 9. We first note that 



1 < X < xm = 2 



1 r-, ^ 

n + - + y^n{n + 1) 



(C16) 



The lower bound comes from r = |,^| > 0, and the upper bound from the requirement that 
the eigenvalues of (which are expectation values) are positive, that is > = sh^(r). 
For the values of 6c such that the discriminant C ^ 0; R{x) of Eq. (|C13b ) is negative for 
X in the interval [x-, x^], where the roots of the polynome are 

^±(0,^) = ^|^ + n5±[(n5)2-|c|%in2(2e)]'^'| . (C17) 
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For a given value of 5, the length of the interval is a strictly decreasing function of 0(^'. 
the smaller the deviation between the super-fluctuant directions of \{v, and p, the 
smaller the allowed range of x. Indeed, the smallest root decreases from x-{0;6) = 
2[n + l/2 + |c|((5)]/(l + 4n(5) > xm to x_(6'm;5) = {l + 2n6)/2A{9M), while the largest root 
decreases from x+(0; 5) = 4[n + 1/2 + |c|(5)] > xm to x+{9m] Sr) = x-{9m'-, Sr) < xm- 

The common value of the roots at 9m means that for this angle, there is a unique value 
of the squeezing parameter. Not surprisingly for 5 = 1, this common value is 1, i.e. r = 0, 
and the fact that 9 can take any value corresponds to the isotropy of the coherent states. 

To study the limit 5 — > 0, please notice that we can set = and that x+(0, 6) is larger 
that Xm for ^ 1. It means that x bellongs to the interval [x-{9),xm] which shrinks to 
Xm since x-{d]5) approaches xm from below. In the limit of a pure state (5 = 0), there 
is a unique representation of the form (|C4p . the state itself (that is the distribution of Eq. 

(IClisPc = 41c(^i))- 
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